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We study the homogeneous and the spatially periodic instabilities in a nematic liquid crystal layer 
subjected to steady plane Couette or Poiseuille flow. The initial director orientation is perpendicular 
to the flow plane. Weak anchoring at the confining plates and the influence of the external electric 
and/or magnetic field are taken into account. Approximate expressions for the critical shear rate 
are presented and compared with semi-analytical solutions in case of Couette flow and numerical 
solutions of the full set of nematodynamic equations for Poiseuille flow. In particular the dependence 
of the type of instability and the threshold on the azimuthal and the polar anchoring strength and 
external fields is analysed. 



I. INTRODUCTION 

Nematic liquid crystals (nematics) represent the sim- 
plest anisotropic fluid. The description of the dynamic 
behavior of the nematics is based on well established 
equations. The description is valid for low molecular 
weight materials as well as nematic polymers. 

The coupling between the preferred molecular orienta- 
tion (director n) and the velocity field leads to interesting 
flow phenomena. The orientational dynamics of nemat- 
ics in flow strongly depends on the sign of the ratio of 
the Leslie viscosity coefficients A = 0,3)0.1. 

In typical low molecular weight nematics A is positive 
(flow- aligning materials). The case of the initial direc- 
tor orientation perpendicular to the flow plane has been 
clarified in classical experiments by Pieranski and Guyon 
Q, Q] and theoretical works of Dubois- Violette and Man- 
neville (for an overview see |3|)- An additional external 
magnetic field could be applied along the initial director 
orientation. In Couette flow and low magnetic field there 
is a homogeneous instability |^{. For high magnetic field 
the homogeneous instability is replaced by a spatially 
periodic one leading to rolls [J. In Poiseuille flow, as 



in Couette flow, the homogeneous instability is replaced 
by a spatially periodic one with increasing magnetic field 
All these instabilities are stationary. 

Some nematics (in particular near a nematic-smectic 
transition) have negative A (non- flow- aligning materials). 
For these materials in steady flow and in the geometry 
where the initial director orientation is perpendicular to 
the flow plane only spatially periodic instabilities are ex- 
pected |5|. These materials demonstrate also tumbling 
motion 6] in the geometry where the initial director 
orientation is perpendicular to the confined plates that 
make the orientational behavior quite complicated. 

Most previous theoretical investigations of the orien- 
tational dynamics of nematics in shear flow were carried 
out under the assumption of strong anchoring of the ne- 
matic molecules at the confining plates. However, it is 
known that there is substantial influence of the bound- 
ary conditions on the dynamical properties of nematics 
in hydrodynamic flow pf la 111 LUj ■ Indeed, the anchoring 
strength strongly influences the orientational behavior 
and dynamic response of nematics under external elec- 
tric and magnetic fields. This changes, for example, the 
switching times in bistable nematic cells 0], which play 
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an important role in applications Recently the in- 
fluence of the boundary anchoring on the homogeneous 
instabilities in steady flow was investigated theoretically 

□ 

In this paper we study the combined action of steady 
flow (Couette and Poiseuille) and external fields (electric 
and magnetic) on the oricntational instabilities of the ne- 
matics with initial orientation perpendicular to the flow 
plane. We focus on flow-aligning nematics. The external 
electric field is applied across the nematic layer and the 
external magnetic field is applied perpendicular to the 
flow plane. We analyse the influence of weak azimuthal 
and polar anchoring and of external fields on both homo- 
geneous and spatially periodic instabilities. 

In section II the formulation of the problem based on 
the standard set of Ericksen-Leslie hydrodynamic equa- 
tions |12| is presented. Boundary conditions and the 
critical Freedericksz field in case of weak anchoring are 
discussed. In section III equations for the homogeneous 
instabilities are presented. Rigorous semi-analytical ex- 
pressions for the critical shear rate a 2 , for Couette flow 
(section III A), the numerical scheme for finding a 2 for 
Poiseuille flow (section III B) and approximate analyt- 
ical expressions for both types of flows (section III C) 
are presented. In section IV the analysis of the spatially 
periodic instabilities is given and in section V we discuss 
the results. In particular we will be interested in the 
boundaries in parameter space (anchoring strengths, ex- 
ternal fields) for the occurrence of the different types of 
instabilities. 

II. BASIC EQUATIONS 

Consider a nematic layer of thickness d sandwiched 
between two infinite parallel plates that provide weak 
anchoring (Fig. ^a). The origin of the Cartesian coordi- 
nates is placed in the middle of the layer with the z axis 
perpendicular to the confining plates (z = ±d/2 for the 
upper/lower plate). The flow is applied along x. Steady 
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FIG. 1: Geometry of NLC cell (a). Couette (b) and Poiseuille 
(c) flows. 

Couette flow is induced by moving the upper plate with 
a constant speed (Fig. \I\b). Steady Poiseuille flow is in- 
duced by applying a constant pressure difference along x 
(Fig. He). An external electric field Eq is applied along 
z and a magnetic field Hq along y. 

The nematodynamic equations have the following form 

□ 

p(d t + v • VM = -p,i + [TV + Tjt] j , (1) 
7l N = -(l-nn-)(72^-n + h), (2) 

where p is the density of the NLC and p , = AP/Ax the 
pressure gradient; 71 = 0:3 — 012 and 72 = «3 + c*2 are 
rotational viscosities; N = n t + v ■ Vn — (V x v) x n/2 
and Aij = (yi j + Vj^)/2, hi = SF/dni. The notation 
/.i = dif is used throughout. The viscous tensor TV- and 
elastic tensor T? are 

TV = ain^jAkmnkUm + a 2 n l N, + a 3 njNi 



J-ij — 



+a 4 A i: j + a 5 nin k A ki + a 6 A lk n k n 3 , (3) 
dF 

(4) 



OUk,i 

where at are the Leslie viscosity coefficients. The bulk 
free energy density F is 

F= i{ifii(V-n) 2 +X 22 [n-(Vxn)] 2 

+%[nx(Vxn)f-£ £ B (n-E„) 2 (5) 
-fi Xa(n ■ H ) 2 j. 
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Here Ku are the elastic constants, e a the anisotropy of 
the dielectric permittivity and \a is the anisotropy of the 
magnetic susceptibility. 

In addition one has the normalization equation 



n = 1 

and incomprcssibility condition 
V-v = 0. 



(6) 



(7) 



The basic state solution of equations (IJ and (0) has the 
following form 



n = (0, 1, 0), v = (v Qx (z), 0, 0), 



(8) 



where vq x = Vq(1/2 + z/d) for Couette and vq x — 
(AP/Ax)[d 2 /a 4 ][l/4 - z 2 /d 2 } for Poiscuillc flow. 

In order to investigate the stability of the basic state 
(JSJl with respect to small perturbations we write: 



n = n + n 1 (z)e^e i(!y , v = v + vi(z)e <Tt e 



<rt iqy. 



(9) 



We do not expect spatial variation along x for steady 
flow. The case q — corresponds to a homogeneous 
instability. Here we analyse stationary bifurcations, thus 
the threshold condition is a = 0. 

Introducing the dimensionless quantities in terms of 
layer thickness d (typical length) and director relaxation 
time t c i = {— a 2 )d 2 / K 22 (typical time) the linearised 
equations and J5J) can be rewritten in the form 

(»7i3 - l)q 2 Sn lz + iqimsq 2 ~ d 2 z )v lx = 0, (10a) 
d z [m2q 2 + (I - V32)d 2 ](Sn lx ) 

+ (Vi2q 4 - Vi2q 2 d 2 z + d 4 z )v ly = 0, (10b) 

{d 2 - k 3 2q 2 - h)ni x + Sn lz + iqv lx = 0, (10c) 
d z (k\ 2 d 2 z - k Z2 q 2 -h + ki 2 e)n lz 

+ \d z (Sn lx ) - {q 2 + Xd 2 )v ly = 0, (lOd) 

vi Z:Z = -iqviy. (10c) 

where = rii/Vj, Vi = («4 + "5 - a 2 )/2, rj 2 = 
(a 3 + a A + a 6 )/2, 773 = a 4 /2, 774 = ct\ + r\\ + T) 2 , tj 5 = 
-(a 2 +a 5 )/2, hj = Ku/Kjj, A = a 3 /a 2 , h = n 2 H$/H%,, 



e = sgn(e a )ir 2 E 2 /E 2 F and H F = (ir/dj^/K^/ifioXa), 
Ep = (tt / d) ^ Ku / (eo\e a \) are the critical Freedericksz 
transition fields for strong anchoring. 

For the shear rate S one has, for Couette flow, 



2 2 VoTd 



b = a , a — 



and for Poiseuille flow 



a 2 2 

S = —a z, a 



Ax 773 



(11) 



(12) 



The anchoring properties are characterised by a surface 
energy per unit area, F s , which has a minimum when 
the director at the surface is oriented along the easy axis 
(parallel to the y axis in our case). A phenomenological 
expression for the surface energy F s can be written in 
terms of an expansion with respect to (n — no). For 
small director deviations from the easy axis one obtains 

F s = \w a n 2 x + \w p n\ z , W a > 0, W p > 0, (13) 

where W a and W p are the "azimuthal" and "polar" an- 
choring strengths, respectively. W a characterizes the sur- 
face energy increase due to distortions within the surface 
plate and W p relates to distortions out of the substrate 
plane. 

The boundary conditions for the director perturbations 
can be obtained from the torques balance equation 



±_ aj, _ + P. = o, 



(14) 



~ d(dnu/dz) driu 
with "±" for z = ±d/2. The boundary conditions Ijl3|l 
can be rewritten in dimensionless form as: 

iftiii/ + n lx = 0, ±f3 p n 1ZtZ + n lz = 0, (15) 

with "±" for z = ±1/2. Here we introduced dimension- 
less anchoring strengths as ratios of the characteristic 
anchoring length (Ka/Wi) over the layer thickness d: 

Pa = K 22 /(W a d), P - K n /(W p d). (16) 

In the limit of strong anchoring, (/3 a , (3 P ) — > 0, one has 
n\ x = n iz = at z = ±1/2. For torque-free boundary 
conditions, (/3 a , (3 P ) — > 00, one has ni XlZ = ni ZyZ = at 
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TABLE I: Symmetry properties of the solutions of equations 
under {z — ► —z}. 



has the form 





Couette flow 


Poiseuille flow 


Perturbation 


"odd" 


"even" 


"odd" 


"even t 


ni x 


odd 


even 


odd 


even 


niz 


odd 


even 


even 


odd 


Vlx 


odd 


even 


odd 


even 


Vly 


even 


odd 


odd 


even 


Viz 


odd 


even 


even 


odd 



the boundaries. From i|16|) one can see that by changing 
the thickness d, the dimensionless anchoring strengths 
(3 a and (3 P can be varied with the ratio f3 a /Pp remaining 
constant. 

The boundary conditions for the velocity field (no-slip) 
are 

= ±1/2) = 0; (17) 
v ly (z = ±1/2) = 0; (18) 
v lz (z = ±1/2) = v lz , z (z = ±1/2) = 0. (19) 

The existence of a nontrivial solution of the linear ordi- 
nary differential equations (|10fl with the boundary condi- 
tions (|15fl . i|17l - H9[l gives values for the shear rate So(q) 
(neutral curve). The critical value S c (q c ), above which 
the basic state (|HJ becomes unstable, are given by the 
minimum of So with respect to q. 

The symmetry properties of the solutions of equations 
(|10|l under the reflection z — > — z is shown in the Table [I] 
We will always classify the solutions by the z symmetry 
of the x component of the director perturbation n\ x (first 
row in Table I). 

In case of positive e a , for some critical value of the 
electric field the basic state loses its stability already in 
the absence of flow (Freedericksz transition). Clearly the 
Freedericksz transition field depends on the polar anchor- 
ing strength. There is competition of the elastic torque 
Kuniz^zz and the field-induced torque SaSoE^niz. The 



n\ z = C cos(irSz/d), 



(20) 



where S = Ef eak /E F and Ef eak is the actual 
Freedericksz field. 

After substituting n\ z into the boundary conditions 
(|15|l we obtain the expression for 5: 



ttS _ 1 

tan T"^- 



(21) 



One easily sees that S — > 1 for (3 p — > and 5 — » 
for f3 p -> 00. For (3 P = 1 one gets Ef eak = 

0A2E F . 

III. HOMOGENEOUS INSTABILITY 

In order to obtain simpler equations we use the renor- 
malized variables as in Ref. 10]: 

S = r x S, N lx = /TV*, N u = n u , V lx = p-'vix, 



Viy = {Fm^viy, Viz = (/Prta)' 1 *!* 
with 



a. 



= "32fc21??32, = — 



Q , 



(22) 



(23) 



In the case of homogeneous perturbations (q = 0) Eqs. 
(fTTJl reduce to V\ z = and 



v ly ,zz-{i-mz){SNi x ), z = Q, 

SN lz - N lx .zz + hN lx = 0, 



(24a) 
(24b) 



solution of Eq. IjlOdfl with n\ x = 0, v\ y = for h 







rfaSNi x + Niz, zz - V ly , z - (k 21 h - e)Ni z = 0. (24c) 

A. Couette flow 

For Couette flow we can obtain the solution of l|24|) 
semi- analytically. For the "odd" solution one gets 

N lx = Ci sinh(az) + C 2 sin(6z), (25) 

N lz = C 3 sinh(C^) + C 4 sinfez), (26) 

Vi y = C 5 cosh(Crz) + C 6 cos(e 2 ^) + C 7 . (27) 
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Taking into account the boundary conditions l|15l 118(1 
the solvability condition for the Ct ( "boundary determi- 
nant" equal to zero) gives an expression for the critical 
shear rate a r : 



For the "even" solution one obtains: 



(/i + d)KiA,cosh(£i/2) 



*(6/2)] 



where 



x [6/3 P cos(6/2)+sin(6/2)] 
(/i-^)[6/3aCos(6/2)+sin(6/2)] 
x [a/3 P cosh(£i/2) + sinhfo/2)] = 0. (28) 



,2 = [(l + fci 2 )ft-fci 2 e]+^ 
41 2fc 12 
, 2 -[(1 + fe 12 )/i - k 12 e] + £ 



2A: 
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,/[(! -fc 12 )fe-fc 12 e] 2 + 4fc 1 2 2 o 4 



(29) 
(30) 
(31) 



JV la = Ci coshfoz) + C* 2 cos(6z) + C 3j (32) 
N lz = C 4 cosh(&*) + C 5 cos(£ 2 z) + C 6 , (33) 
V ly = C 7 smh(Ciz) + C 8 z. (34) 



The boundary conditions 1(17119(1 now lead to the fol- 
lowing condition ("boundary determinant"): 



-6/3aSin(6/2) + cos(6/2) (h + sin(£ 2 /2) + cosfe/2)] 

Ci/3 a sinh(a/2) + cosh(&/2) (ft - sinli(£i/2) + cosh(fc/2)] 



7y 23 ^ - k 12 e) _ 
a 4 fci2?723 
• 23 sin(£ 2 /2) 



1 6 

1 - V23 



sinh(Ci/2) 



(35) 



r 



B. Poiseuille flow 

In the case of Poiseuille flow the system ((24(1 with S = 
—za 2 //3 admits an analytical solution only in the absence 
of external fields (in terms of Airy functions) [lfj . In the 
presence of fields we solve the problem numerically. In 
the framework of the Galerkin method we expand Ni x , 
Ni z and V\ y in a series 



N lx = C ltn f n (z), 

71=1 
OO 

n=l 

oo 

Vly = ^ C 3,n«n(z), 



(36) 



where the trial functions /„ , g n and u n satisfy the bound- 
ary conditions ((15|) , ((18(1 . For the "odd" solution we write 



fn{z) = Q(z;P a ), g n (z) = Q(z;0 p ), u n (z) = v° n {z) 



(37) 



and for the "even" solution 



fn(z) = Q(z;fi a ), g n (z) = C(z;P p ), u n (z) = v e n (z). 

(38) 

The functions £°(z;/3), £®(#;/3), Vn(z), v^{z) are given in 
Appendix A. In our calculations we have to truncate the 
expansions ((36(1 to a finite number of modes. 

After substituting 1(36(1 into the system 1(24(1 and pro- 
jecting the equations on the trial functions f n {z), g n {z) 
and u n [z) one gets a system of linear homogeneous alge- 
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braic equations for X = {Ci )H } in the form (A— a 2 £>)X = 
0. We have solved this eigenvalue problem for a 2 . The 
lowest (real) eigenvalue corresponds to the critical shear 
rate a 2 . According to the two types of z-symmetry of 
the solutions (and of the set of trial functions) one ob- 
tains the threshold values of a 2 , for the "odd" and "even" 
instability modes. The number of Galerkin modes was 
chosen such that the accuracy of the calculated eigenval- 
ues was better than 1% (wc took ten modes in case of 
"odd" solution and five modes for "even" solution). 

C. Approximate analytical expression for the 
critical shear rate 

In order to obtain an easy-to-use analytical expres- 
sion for the critical shear rate as a function of the sur- 
face anchoring strengths and the external fields we use 
the lowest-mode approximation in the framework of the 
Galerkin method. By integrating (|24a[) over z one can 
eliminate V\ v , z from Q24c|l which gives 

SN lx + N 1ZiZZ + {k 2 ih-e)N lz = K, (39) 

where K is an integration constant. Taking into account 
the boundary conditions for V\ y one has 

1/2 

K-(l- m2 ) J SN lx (z)dz = 0. (40) 

-1/2 

We choose for the director components Ni x , N% z the 
one-mode approximation 

N lx = df(z), N lz = C 2 g(z), (41) 

Substituting (|4*TJl into (|24b|l and JSHJ and projecting the 
first equation on f(z) and the second one on g(z) we 
get algebraic equations for C{. The solvability condition 
[together with l|40[l ] gives the expression for the critical 
shear rate a 2 

«l - ^ (42) 



TABLE II: Trial functions for the homogeneous solutions. 





Couette flow Poiseuille flow 


Function 


"odd" "even" "odd" 


"even" 


/(*) 


C?(z; p a ) a(z; M Cf(»; AO 


Cf(*; Pa) 




Cf(*; P v ) Pv) Ci e («; Pv) 


<»1 ( z > Pp) 


where c\ 


= {ff")-Hf),c 2 = (gg")-(h/ku 


-e){g 2 ), 


ca = (sfg)[{sfg) - (1 - %3) (sf)(g)], where (. . . 


) denotes 


a spatial 


average 






1/2 






(...)= J {...)dz. 


(43) 



-1/2 



The values for the integrals (...) are given in Appendix 
B. In Table ITT1 and Appendix A the trial functions used 
are given. Equation (|42|l can be used for both Couette 
and Poiseuille flow by choosing the function s(z) [where 
s(z) = 1 for Couette flow and s(z) = —z for Poiseuille 
flow] and the trial functions f(z) and g(z) with appro- 
priate symmetry. 

For the material MBBA in the case of Couette flow 
the one-mode approximation l|42|l for the "odd" solution 
gives an error that varies from 2.5% to 16% when Ho/Hp 
varies from to 4. The "even" solution has an error of 
0.6% 8% for sC Hq/H f ^ 3 and of 0.6% 12% for 
< E /E F < 0.6. 

For Poiseuille flow for "odd" solution the error is 29% 
in the absence of fields. For the "even" solution the error 
is 12% ~ 15% for magnetic fields < H /H F ^ 0.5. 

For both Couette and Poiseuille flow the accuracy of 
the formula (|42J) decreases with increasing field strengths. 

IV. SPATIALLY PERIODIC INSTABILITIES 

We used for Eqs. (|10fl again the renormalized variables 
Q22JI. The system (I10|) has no analytical solution. Thus 
we solved the problem numerically in the framework of 
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TABLE III: Trial functions for the spatially periodic solutions. 





Couette flow 


Poiseui 


He flow 


Him c{~ i nn 


"odd" 


"even" 


"odd" 




/(*) 








£(*; Pa) 


<?(*) 


ft) 


ft,) 


C5(z; ft.) 


ft) 


u{z) 










w(z) 


«£(*) 








the Galerkin method: 










oo 




OO 




N lx = e*« Y Ci, n fn(z), N lz = 


= e 1 ™ Y °2,7 


l5 „(z),(44) 




71 = 1 




n=l 






OO 




OO 




V lx = Y C 3 ,nu n 


(«), ^ - 


■■e^Y C ^ n {z){^>) 




n=l 




71=1 




After 


substituting 


gg) 


into the 


system 


TO and 


projecting 


on to 


the trial 


functions 


{/n(z), 5n 


(z), u„(z), 




we get a 


system of 



linear homogeneous algebraic equations for X = {Ci >n }. 
This system has the form [A(q) - a 2 (q)B(q)]'X — 0. We 
have solved the eigenvalue problem numerically to find 
the marginal stability curve a{q). For the numerical 
calculations we have chosen the trial functions shown in 
Table |inl and Appendix A. 

In order to get an approximate expression for the 
threshold we use the leading-mode approximation in the 
framework of the Galerkin method. We used the same 
scheme described above for the single mode and get the 
following formula for the critical shear rate: 

al = \Zv23fxfz/(a 2 a3), (46) 

with 

f x = (ff")-(q 2 k 3 2 + h)(f), (47) 

U = (.9.9") - (q 2 hi + k 12 h - e)(g 2 ), (48) 

&2 = [{fag) - q 2 (l - mi)(fu)(gsu)fr], (49) 

"3 = (fsg) + [a 23 q 2 (gw) + a 3 (gw")} (50) 

x [(1 - m2)(w[sf}") - V52 q 2 (wsf)]/r, (51) 

1 = q 2 {uu) -mi{uu"), (52) 

r = (ww^) — r\viq 2 '(ww") + rj 12 q 4 (ww) . (53) 



The values of the integrals (...) appearing in the expres- 
sion 146(1 are given in Appendix C. 

In the case of strong anchoring an approximate analyti- 
cal expression for a 2 = a 2 (q c ) was obtained by Manneville 
using test functions that satisfy free-slip boundary 
conditions. The formula H4(jfl is more accurate because 
we chose for v\ z Chandrasekhar functions that satisfy the 
boundary conditions 1(19(1 . 

For calculations we used material parameters for 
MBBA. The accuracy of J2BJ is better than 1% for Cou- 
ette flow and better than 3% for Poiseuillc flow. Note, 
that Eq. I(42() for the homogeneous instability is more 
accurate than 1(46(1 for q = because 1(46(1 was obtained 
by solving four equations 1(11)1) by approximating all vari- 
ables, whereas ((4211 was obtained by solving the reduced 
equations ((24(1 by approximating only two variables. 

V. DISCUSSION 

For the calculations we used parameters for MBBA 
at 25 °C Calculations were made for the range of 

anchoring strengths ft, = -j- 1 and ft, = -j- 1 . 

A. Couette flow 

We found that without and with an additional electric 
field the critical shear rate a 2 for the "even" type ho- 
mogeneous instability (EH) is systematically lower than 
the threshold for other types of instability (Fig. |5Ji-c). 
Note, that in the presence of the field the symmetry with 
respect to the exchange ft, <-> ft, is broken. 

In Fig. [3 contour plots for the critical value a 2 vs. 
anchoring strengths ft, and ft, for different values of the 
electric field are shown. The differences between a 2 ob- 
tained from the exact, semi-analytical solution 1(35(1 and 
from the one-mode approximation ((42|) are indistinguish- 
able in the figure. 

In Fig. contour plots of a 2 (thin dashed lines) and 
the boundaries where the type of instability changes [the 
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FIG. 2: Contour plot of the critical shear rate 0% for Couette 
flow vs. p a and P P . a: E = 0; b: E = E% eak , e a < 0; 
c: E = Ep eak , e a > 0. E F is defined after Eq. J2J and 
calculated in Eq. H2H . 



FIG. 3: Critical shear rates and phase diagram for the insta- 
bilities under Couette flow with additional magnetic field, a: 
H /Hf = 3; b: Ho/ Hf = 3.5; c: H /H F = 4. Boundaries for 
occurrence of instabilities are given by thick solid lines (full 
numerical) and thick dashed lines (one-mode approximation) . 
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solid lines are obtained numerically, the thick dashed 
lines from Q42[l] for different values of magnetic field are 
shown. For not too strong magnetic field in the region 
of weak anchoring the "odd" type homogeneous insta- 
bility (OH) takes place (Fig. Eli). Increasing the mag- 
netic field the OH region expands toward stronger an- 
choring strengths. Above Hq rs 3.2 a region with lowest 
threshold corresponding to the "even" roll mode (ER) 
appears. This region has borders with both types of the 
homogeneous instability (Fig. Ob). With increasing mag- 
netic field the ER region increases (Fig. [jj) and above 
Hq/Hf — 4 the ER instability has invaded the whole 
investigated parameter range. For strong anchoring and 
Hq/Hf = 3.5 the critical wave vector is q c — 5.5. It in- 
creases with increasing magnetic field and decreases with 
decreasing anchoring strengths. With increasing mag- 
netic field the threshold for the EH instability becomes 
less sensitive to the surface anchoring. Leslie has pointed 
out (using an approximate analytical approach) that for 
strong anchoring a transition from a homogeneous state 
without transverse flow (EH) to one with such flow (OH) 
as the magnetic field is increased is not possible in MBBA 
because of the appearance of the ER type instability • 
This is consistent with our results. We find that the EH- 
OH transition in MBBA is possible only in the region of 
weak anchoring (Figs. |3Ji,-c). 

In Fig. 0] marginal stability curves for different values 
of the magnetic field and fixed anchoring strengths is 
shown (solid line for ER and dashed lines for OR). There 
are always two minima for the even mode; one of them 
at q = that corresponds to the homogeneous instability 
EH. For small magnetic field the absolute minimum is 
at q = (line a). The OR curve is systematically higher 
than ER. In a small range of q (dotted lines) a stationary 
ER solution does not exist but we have OR instead. With 
increasing magnetic field the critical amplitude for the 
EH minimum (q = 0) increases more rapidly then the one 
for the ER minimum (q ^ 0) so that for Hq/Hf > 3.4 
the ER solution is realized (lines b and c). The range of 
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FIG. 4: a c vs. q. Couette flow, /3 a — 0.1, f3 p = 0.1. a: 
H / H F = 3; b: Hq/Hf = 3.4; c: H /H F = 4. 

q where ER is replaced by OR expands with increasing 
magnetic field. 

For the ER instability in the absence of fields and 
strong anchoring we find a 2 c — 12.15 from the semi- 
analytical expression l|35|) as well as from the one-mode 
approximation l|42l) and also (|46|l with q = 0. The only 
available experimental value for a 2 , is 6.3 ± 0.3 [2|. We 
suspect that the discrepancy is due to deviations from the 
strong anchoring limit and the difference in the material 
parameters of the substance used in the experiment. As- 
suming [3 a <C 1 one would need (3 P w 1 to explain the 
experimental value. 

B. Poiseuille flow 

In Fig. [S] the contour plot for a 2 c [thin dashed lines 
from the full numerical calculation, dotted lines from the 
one-mode approximations (|42|l and (|46|l ] and the bound- 
ary for the various types of instabilities [thick solid line: 
numerical; thick dashed line: 142(1 and l|46|) ] are shown. 
In Poiseuille flow the phase diagram is already very rich 
in the absence of external fields. In the region of large 
j3 a one has the EH instability. For intermediate anchor- 
ing strengths rolls of type OR occur [Fig. [St], Note, 
that even in the absence of the field there is no symme- 
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FIG. 5: Critical shear rates and phase diagram for the in- 
stabilities in Poiseuille flow, a: E = 0; b: E = E^ ak , 
e a < 0; c: Eo = Eo eah , e a > 0. Thin dashed lines: full 
numerical threshold; dotted lines: one-mode approximation 
for threshold. Boundaries for occurrence of instabilities are 
given by thick solid lines (full numerical) and thick dashed 
lines (one-mode approximation). 
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FIG. 6: Phase diagram for the instabilities under Poiseuille 
flow with an additional magnetic field (Ho/Hf = 0.4). 

try under exchange (3 a P , contrary to Couette flow. 
The one-mode approximations 142J1 and 146(1 not give the 
transition to EH for strong anchoring. Here we should 
note that in that region the difference between the EH 
and the OR instability thresholds is only about 5%. By 
varying material parameters [increase ai by 10% or de- 
crease q;3 by 20% or as by 25% or K 33 by 35%] it is 
possible to change the type of instability in that region. 

Application of an electric field leads for e a < (e a > 0) 
to expansion (contraction) of the EH region [Figs. |5j> and 
Et]. At Eq/Ef = 1 and e a < rolls vanish completely 
and the EH instability occurs in the whole area investi- 
gated. For e a > the instability of OH type appears in 
the region of large (3 p . In this case, increasing the electric 
field from Ep eak to Ep cause an expansion of the OH re- 
gion. Note that for (3 P > 1, which is in the OH region, 
the Freedericksz transition occurs first . 

An additional magnetic field suppresses the homoge- 
neous instability (Fig. EJl. Above H /Hp w 0.5 the OR 
instability (Fig. occurs for all anchoring strengths in- 
vestigated. 

The wave vector q c in the absence of fields is 1.4. Ap- 
plication of an electric field decreases q c whereas the mag- 
netic field increases q c . The wave vector decreases with 
decreasing anchoring strengths. 
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In the absence of fields and strong anchoring we find 
for the EH instability a c = 102 [Eq. (JE) gives 110 and 
Eg. l|46l) with q — gives 130]. The experimental value is 
92 |l6J. Thus, theoretical calculations and experimental 
results are in good agreement. Note, that in the exper- 
iments 16] actually not steady but oscillatory flow with 
very low frequency was used (/ = 5 ■ 10~ 3 Hz). 

In summary, the orientational instabilities for both 
steady Couette (semi-analytical for homogeneous insta- 
bility and numerical for rolls) and Poiseuille flow (numer- 
ical) were analysed rigorously taking into account weak 
anchoring and the influence of external fields. Easy-to- 
use expressions for the threshold of all possible types of 
instabilities were obtained and compared with the rigor- 
ous calculations. In particular the region in parameter 
space where the different types of instabilities occurred 
were determined . 
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APPENDIX A: TRIAL FUNCTIONS 



APPENDIX B: INTEGRALS FOR THE 
HOMOGENEOUS INSTABILITY 

1. Couette flow 

"Odd" solution: (*/) = (g) = 0, (f 2 ) = (3 + S2(3 a + 
12^ 2 /3 2 )/6, (g 2 ) = (3 + 32/3 p + 12^ 2 /3 2 )/6, (sfg) = 
[3 + 16(/3„ + f3 p ) + 127r 2 /3 a /3 p ]/6, (//") = -2(3 + 20{3 a + 
3tt 2 /52)/3, (gg") = -2(3 + 20/3 p + 3^ 2 /3 2 )/3. 
"Even" solution: (sf) = (2 + tt 2 f^ a ) / n , (g) = (2 + 
7r 2 (3 p )/n, </ 2 ) = (1 + 8/3 a + 2tt 2 /3 2 )/2, (g 2 ) = (1 + 
8f3 p + 2ir 2 /3 2 )/2, (sfg) = [1 + 4(/3 a + (3 p ) + 2tt 2 /? q /3 p ]/2, 
(//") = tt 2 (1 + 4fl,)/2, (gg") = tt 2 (1 + 4/3 p )/2. 

2. Poiseuille flow 

"Odd" solution: (sf) = -(1 + 8/? q )/(2tt), (g) = 
-(2 + 7t 2 /3 p )/tt, (f 2 ) = (3 + 32[3 a + 12^ 2 /3 2 )/6, (g 2 ) = 
(1 + 8/3 p + 2tt 2 /3 2 )/2, (sfg) = -[16 + 97r 2 (/? a + /3 P ) + 
727r 2 /3 a /5 p ]/(187r 2 ), (//") = -27r 2 (3 + 20/3 a + 3^ 2 /3 2 )/3, 
(gg") = -^(l + A0 p )/2. 

"Even" solution: (sf) = 1(g) = 0, (f 2 ) = (1 + 8/3 Q + 
2tt 2 /3 2 )/2, (g 2 ) - (3 + 32f3 p + 12tt 2 /3 2 )/6, (sfg) = 
-[16+97r 2 (/3 Q +/3 p )+72^ 2 /3 a /3 p ]/(187r 2 ), (//") = -tt 2 (1 + 
4/3 a )/2, (gg") = -2^ 2 (3 + 20/3 p + 3tt 2 /? 2 ). 



In the calculations we used the following set of trial 
functions: 

Q(z;j3) = sin(2mrz) + 2nir/3 sin([2n - l]irz), 

Q(z;P) = cos([2n - l]nz) + (2n - 1)tt/3 cos(2[n - l]wz), 

v° L (z) = sm(2nnz), V^{z) — cos([2n — 1]ttz), 

sinh(A 2 „z) _ sin(A 2 »2:) 

in[Z) ~ sinh(A 2 „/2) sin(A 2n /2)' 

e cosh(A 2 „_iz) _ cos(A 2 „_iz) 

^ n{Z > ~ cosh(A 2n -i/2) cos(A 2 „_i/2) ' 

$n(z) and <;!„(z) are the Chandrasekhar functions and A„ 
are the roots of the corresponding characteristic equa- 
tions |l7j . 



APPENDIX C: INTEGRALS FOR THE 
SPATIALLY PERIODIC INSTABILITY 

1. Couette flow 

"Odd" solution: (wsf) « 0.69043 + 3.2870/3 Q , (w[sf]") « 
-27.258 - 32.441/3 a , (/ 2 ) = (3 + 32/3 a + 12tt 2 /3 2 )/6, 
(//") - -^ 2 (6 + 40/3 Q + 6tt 2 /3 2 )/3, (/*<?) - (3 + 16(/? + 
/3 P ) + 127T 2 /3 a /3 p )/6, ( 5 su) = (3 + 16/3 p )/6, ( 5 2 ) = (3 + 
32/3 p + 12tt 2 /5 2 )/6, < 55 ") = -tt 2 (6 + 40/3 p + 6tt 2 /3 2 )/3, 
(u 2 ) = 1/2, (W) = -2vr 2 , (fu) = (3 + 16/3 a )/6, 
(w 2 ) = 1, (ww") w -46.050, (W 4 )) = 3803.5, (gw) « 
0.69043 + 3. 2870/3 p , (gw") w -27.257 - 32.441/3 P . 
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"Even" solution: (wsf) w 0.69739+2.6102/? a , (Ms/]") « 
-6.8828, (/ 2 ) = (1 + 8/3 a + 2tt 2 /3 2 )/2, (//") = -tt 2 (1 + 
4/?„)/2, (,/sg) = (1 + 4(/3 a + /3 P ) + 2^ 2 /3 a /3 p )/2, foau) = 
(l + 4/? p )/2, ( 5 2 ) = (1 + 8/3 p + 2tt 2 /3 2 )/2, (gg") = -tt 2 (1 + 
4/3 p )/2, ( u 2 ) = 1/2, (««") = -2tt 2 , (fu) = (1 + 4/3 a )/2, 
(w 2 ) = 1, (W) w -12.303, (ww^) w 500.56, (gw) w 
0.69738 + 2.6102/? p , (#«/'} « -6.8828. 

2. Poiseuille flow 

"Odd" solution: (wsf) w -0.10292 - 0.49816/3 a , 
<«;[*/]"> « -0.87673 - 22.615/3 a , (/ 2 ) = (3 + 32/3 Q + 
12^ 2 /? 2 )/6, (//") = -7r 2 (6 + 40/? a + 6^ 2 /3 2 )/3, (fag) = 



-(16 + 9^ 2 (/? a +/3 p ) + 727r 2 /3 a /? p )/(18^ 2 ), (gsu) = -(16 + 
9^ 2 /? p )/(18^ 2 ), (g 2 ) = (1 + 8/? p + 2tt 2 /3 2 )/2, (gg") = 
-tt 2 (1 + 4/? p )/2, (u 2 ) = 1/2, (W) = -2tt 2 , (fu) = (3 + 
16/3 a )/6, (w 2 ) = 1, (ww") « -12.303, (W 4 )) w 500.56, 
(.gw) « 0.69738 + 2. 6102/3 p , w -6.8828. 

"Even" solution: (wsf) « -0.12206 - 0.59694/3 a , 
<«;[*/]"> « 4.4917, (/ 2 ) = (1 + 8[3 a + 2tt 2 /3 2 ), (//") = 
-tt 2 (1 + 4/3 a )/2, (/s 5 ) = -(16 + 97r 2 (/? a + /? p ) + 
72^ 2 /3 a /3 p )/(187r 2 ), (gsu) = -(16 + 9^ 2 /? p )/(18^ 2 ), 
(, 9 2 ) = (3 + 32/3 p + 12tt 2 /? 2 )/6, - -2tt 2 (3 + 20/3 p + 

3^ 2 /3 2 )/3, ( U 2 > = 1/2, (uu") = -tt 2 /2, (fu) = (1 + 
4/3 a )/2, (w 2 ) = 1, (ww") w -46.050, (W 4 )} « 3803.5, 
(.gio) « 0.69043 + 3. 2870/3 p , « -27.257- 32.441/3 p . 
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